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ABSTRACT 

Supernova (SN) rates are a potentially powerful diagnostic of star formation history (SFH), metal 
enrichment, and SN physics, particularly in galaxy clusters with their deep, metal-retaining potentials, 
and simple SFH. However, a low-redshift cluster SN rate has never been published. We derive the 
SN rate in galaxy clusters at 0.06 < z < 0.19, based on type la supernovae (SNe la) that were 
discovered by the Wise Observatory Optical Transient Survey. As described in a separate paper, a 
sample of 140 rich Abell clusters was monitored, in which six cluster SNe la were found and confirmed 
spectroscopically. Here, we determine the SN detection efficiencies of the individual survey images, 
and combine the efficiencies with the known spectral properties of SNe la to calculate the effective 
visibility time of the survey. The cluster stellar luminosities are measured from the Sloan Digital Sky 
Survey (SDSS) database in the griz SDSS bands. Uncertainties are estimated using Monte-Carlo 
simulations in which all input parameters are allowed to vary over their known distributions. We 
derive SN rates normalized by stellar luminosity, in SNU units (SNe per century per 10 10 L©) in 
five photometric bandpasses, of 0.36t^4 ± 0.02 (B), 0.351±£^ ± 0.020 (g), 0.288±£;i^ ± 0.018 (r), 
0.229;^ ±0.014 (i), 0.186±£; u7 4± 0.010 (z), where the quoted errors are statistical and systematic, 
respectively. The SN rate per stellar mass unit, derived using a color-luminosity-mass relation, is 
0.098±°°^ ± 0.009 SNe (century IO^Mq)" 1 . The low cluster SN rates we find are similar to, and 
consistent with, the SN la rate in local elliptical galaxies. 

Subject headings: galaxies: clusters: general — supernovae: general — supernovae: individual (SN 
1998eu, SN 1998fc, SN 1999cg, SN 1999ch, SN 1999ci, SN 1999ct) 



1. INTRODUCTION 

The event-rate history of supernovae (SNe), their pro- 
genitor systems, and the physical mechanisms behind 
these explosions, all have important implications for sev- 
eral branches of astrophysics and cosmology. As the 
primary sources of iron and other heavy elements, SNe 
have a key role in the chemical evolution of the Uni- 
verse (e.g., Nomoto et al. 2005). The history of enrich- 
ment is determined by the rate of SNe in the past and 
by the inventory of heavy elements that is released by 
each explosion. In cosmology, type la SNe (SNe la) used 
as cosmological distance indicators have provided direct 
evidence for an accelerated expansion of the Universe 
(e.g., Riess et al. 1998, 2004; Perlmutter et al. 1999; 
Astier et al. 2006). The use of SNe la as standard can- 
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dies relies on the assumptions that there is no evolution 
in their intrinsic properties, and that their luminosities 
can be calibrated through an unevolving empirical re- 
lation to their light-curve shape (for reviews, see Lei- 
bundgut 2001; Filippenko 2004, 2005). If these assump- 
tions could be founded on more solid observational and 
physical grounds, the cosmological use of SNe la would 
stand on a firmer footing. 

It is widely agreed that SNe la occur when a white 
dwarf in a binary system is pushed toward the Chan- 
drasekhar mass limit and undergoes a runaway ther- 
monuclear explosion. However, there is no consensus 
regarding how this situation is reached. In the single- 
degenerate scenario, a white dwarf accretes matter from 
a main-sequence, subgiant, or giant star. In the double- 
degenerate scenario, two white dwarfs merge. Further- 
more, it is not clear how the explosion proceeds. Since 
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direct observational identification of a SN la progenitor 
system is unlikely, indirect evidence may be the only 
practical way to distinguish between models. Modeling 
of the observed spectra of SNe provides many clues (e.g., 
Nomoto et al. 1984), as do studies of Galactic SN rem- 
nants (e.g., Badenes et al. 2006). However, these studies 
mostly constrain the explosion physics (deflagration or 
detonation) and are less sensitive to the identity of the 
mass donor. 

A complementary approach exploits the fact that each 
of the above scenarios may predict a "time delay distribu- 
tion," which is the distribution of times elapsed between 
the formation of a stellar population and the explosions 
of some members of the population as SNe la (e.g., Yun- 
gelson & Livio 2002; Greggio 2005; Forster et al. 2006). 
The time delay, in turn, convolved with the star forma- 
tion history (SFH), dictates the evolution of the SN rate 
with cosmic time or with redshift. By measuring the 
SN la rate vs. redshift z, one can, in principle, constrain 
the form of the time delay distribution and obtain clues 
about the progenitor scenario (e.g., Pain et al. 2002; 
Gal- Yam & Maoz 2004; Strolger et al. 2004; Sullivan et 
al. 2006; Barris & Tonry 2006; Neill et al. 2006; Forster 
et al. 2006). 

In practice, there are several complications. A pre- 
diction of the SN rate is sensitive to the assumed SFH. 
For example, Forster et al. (2006) recently re-analyzed 
the SN la sample of Strolger et al. (2004) and argued 
that, by varying the choice of SFH within the uncer- 
tainties allowed by observations, a wide range of time 
delay distributions are consistent with the SN data, and 
hence none of the SN la scenarios can be ruled out. An- 
other complication follows from the recent evidence that 
SNe la are descended from two distinct stellar popula- 
tions (Mannucci et al. 2005; Scannapieco & Bildsten 
2005; Sullivan et al. 2006), perhaps leading to two SN la 
"channels" - a "prompt" channel, with a delay time of 
order 1 Gyr or less, and a "delayed" channel, with delay 
times that extend to > 10 Gyrs. Effectively, this means a 
time delay distribution that is more complex than those 
predicted by single-channel models. The combination of 
these two issues, combined with debates regarding sam- 
ple completeness, have prevented conclusive determina- 
tions of delay times. 

Measuring the SN la delay time in galaxy clusters, 
however, may provide a solution to these difficulties. 
Analysis of the stellar population of cluster galaxies sug- 
gests that the SFH in clusters is simpler than that in the 
field, and can likely be approximated by a brief starburst 
at high redshift (e.g., Wuyts et al. 2004; Holden et al. 
2005). Sullivan et al. (2006) show that the mix between 
the delayed and prompt components of the SN la pop- 
ulation changes in different environments. The prompt 
component dominates late-type galaxies, while elliptical 
galaxies contain only members of the delayed group. SN 
rate measurements in clusters could then provide an esti- 
mate of the delayed component alone with minimal con- 
tamination. If the population of SNe la is indeed bi- 
modal, pinning down one of its components from cluster 
SN studies could help to disentangle the two populations 
in the field. 

Furthermore, measurement of cluster SN rates may 
help resolve the question of the dominant source of the 
high metallicity in the intracluster medium (e.g., Ren- 



zini et al. 1993; Maoz & Gal- Yam 2004) - SNe la, or 
core-collapse SNe from an early stellar population with 
a top-heavy initial mass function (e.g., Maoz & Gal- Yam 
2004; Loewenstein 2006). Clusters are excellent labora- 
tories for studying enrichment, due to their simple SFH 
and their deep potentials from which matter cannot es- 
cape. Therefore, understanding the source of intracluster 
metals would be relevant for understanding metal enrich- 
ment in general, and the possible role of first-generation 
stars in early-universe enrichment. 

Although SN la rates have been measured extensively 
out to z ~ 1 (Pain et al. 1996, 2002; Cappellaro et al. 
1999; Hardin et al. 2000; Tonry et al. 2003; Madgwick 
et al. 2003; Blanc et al. 2004; Dahlen et al. 2004; Bar- 
ris & Tonry 2006; Neill et al. 2006), there have been 
few attempts to measure the SN rate in galaxy clus- 
ters (e.g., Barbon 1978; Norgaard-Nielsen et al. 1989; 
Reiss et al. 1998). The only published cluster SN rate 
at distances beyond the Virgo cluster is by Gal- Yam, 
Maoz, & Sharon (2002). This measurement was based 
on the detection of two or three likely cluster SNe la in 
archival Hubble Space Telescope (HST) images of high- 
redshift clusters. The derived rates, 0.39±°;25 SNu B 
at z = 0.25 and 0.80±£;;g SNu B at z = 0.9 [1 SNu B 
= 1 SN century _1 (10 10 LB.o) _1 ], have large uncertain- 
ties due to the small number of detected SNe and the 
lack of follow-up observations, a consequence of search- 
ing archival data. 

We are carrying out a program to measure the cluster 
SN rate at both low and high redshifts. In this paper, 
we derive the cluster SN rate based on a low-redshift 
(0.06 < z < 0.19) cluster SN survey, the Wise Observa- 
tory Optical Transient Search (WOOTS). 

The cluster sample, the observational design, and the 
SN detection, spectroscopic follow-up observations, and 
classification are described in a separate paper by Gal- 
Yam et al. (2006), hereafter Paper I. 

Here, we present the details of the derivation of the 
cluster SN rate and briefly discuss its implications. 

Throughout the paper we assume a flat cosmology, 
with parameters J7a — 0.7, J1m = 0.3, and Hq = 
70 km s^Mpc -1 . 

2. THE SURVEY 

In Paper I, we present the observational details of 
WOOTS, the SNe that were discovered by it, and their 
followup observations. Briefly, WOOTS was a survey 
for SNe and other transient or variable sources in the 
fields of 161 galaxy clusters, and was conducted in 1997- 
1999. The cluster sample was selected from the catalog 
of Abell, Corwin & Olowin (1984), based on the follow- 
ing criteria: cluster redshift in the range 0.06 < z < 0.2 
(based on the compilation of Struble & Rood 1991); dec- 
lination S > 0; Abell richness class R > 1 (Abell galaxy 
count N > 65); and cluster radius (estimated by Leir & 
Van den Bergh 1977) smaller than 20". Monthly dark- 
time observations used the Wise Observatory 1 m tele- 
scope, with a 1024 x 1024, 0'.'7 pixel, Tektronics CCD 
imager, giving a 12' x 12' field of view. Images were 
unfiltered (see also § [3|) reaching a limiting magnitude 
equivalent to R ~ 22. Images (of ~ 40 clusters per 
night) were compared with previously obtained template 
images using image subtraction methods, and searched 
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for transients. 

All candidate transients were followed up photomet- 
rically at Wise, and spectroscopically with larger tele- 
scopes. The determination of SN redshifts, ages and 
types was done by comparing the observed spectra to red- 
shifted versions of high S /N template spectra of nearby 
SNe, drawn from the spectroscopic archive presented by 
Poznanski ct al. (2002). In the course of the survey, 
12 SNe were discovered, all of them spectroscopically 
confirmed (11 as type-la, and one as a type-HP SN). 
Seven of the SNe (all of them SNe la) were in their re- 
spective clusters, and the remaining five were foreground 
and background events. All additional transient/ variable 
candidates were confirmed as non-SN events: asteroids, 
active galactic nuclei, and variable stars. 

The SN sample we will analyze in this paper does not 
include SN 2001al (Gal Yam et al. 2003), a cluster SN la 
which was discovered during an extension of WOOTS us- 
ing a new CCD camera during its commissioning phase. 
We do not consider in this paper any of the survey images 
from that survey extension. In addition, we do not in- 
clude in our analysis a sub-sample of 21 clusters for which 
the imaging data were inadvertently lost. Our sample 
thus consists of six confirmed cluster SNe la that were 
discovered by monitoring, effectively, 140 Abell clusters 
at 0.06 < z < 0.19. Table [T] lists the six events and their 
main parameters. 



TABLE 1 
WOOTS Cluster SNe Ia 



SN 


Cluster 


Rcdshift 


1998fc 


Abell 403 


0.10 


1998eu 


Abell 125 


0.18 


1999cg 


Abell 1607 


0.14 


1999ch 


Abell 2235 


0.15 


1999ci 


Abell 1984 


0.12 


1999ct 


Abell 1697 


0.18 



3. PHOTOMETRIC CALIBRATION 

A derivation of the SN rate from our survey requires 
knowledge of the depth of all the survey images, and 
hence a photometric calibration of the survey data in the 
observed bandpass. The WOOTS images were obtained 
in unfiltered ("clear") mode, resulting in a very broad, 
nonstandard photometric bandpass. We have deter- 
mined the effective WOOTS-clear bandpass by multiply- 
ing the quantum efficiency curve of the back-illuminated, 
Lumogen-coated, Tektronix CCD used in the survey, by 
the typical atmospheric transmission at Wise Observa- 
tory for the mean airmass (1.2) of the WOOTS obser- 
vations. We have verified that the effects of varying 
airmass and of mirror reflectivity vs. wavelength have 
negligible effect of the shape of the bandpass, The peak 
transmission of the effective WOOTS-clear band is at 
~ 7000 A, with a full-width at half maximum (FWHM) 
of - 4700 A. 

To obtain the individual photometric zeropoints of all 
the WOOTS images, we begin by using those images 
that fall within the coverage of the Sloan Digital Sky 
Survey Data Release 4 (SDSS DR4; Adelman-McCarthy 
et al. 2006). We extract the SDSS photometry of ob- 
jects identified in each image and find the best fit to the 



SDSS ugriz magnitudes from among a set of template 
spectra of stars (with spectral types O to M; Silva et al. 
1992) and galaxies (Kinney et al. 1996). We then use 
the WOOTS-clear bandpass to perform synthetic, Vega- 
based photometry on the best-fit spectrum for each ob- 
ject. Finally, we compare these synthetic magnitudes to 
the instrumental magnitudes of the same objects, defined 
as 77ij nst = — 2.51og(counts). The WOOTS zeropoint is 
then defined as the median of the differences between 
the synthetic and the instrumental magnitudes of the 
calibration sources in the image. The root-mean-squarc 
(rms) scatter of these differences defines the statistical 
zeropoint error for each image. Every source in an image 
can then be calibrated to a Vega-based WOOTS-clear 
magnitude by adding to its instrumental magnitude the 
derived zeropoint of the image. We estimate the accuracy 
of the zeropoints measured this way by adding the sys- 
tematic uncertainty in the SDSS zeropoints (0.02 mag) 
to the statistical zeropoint error determined for each im- 
age, as defined above (typically also 0.02 mag), resulting 
in a total zeropoint error of 0.04 mag. 

In order to calibrate the images that are in areas not 
covered by the SDSS, we compare the instrumental mag- 
nitudes of objects in our images with their "RV mag- 
nitudes in the USNO-B catalog (Monet et al. 2003), to 
get an "i?l"-based zeropoint. The USNO-i?l band has 
a similar effective wavelength to that of the WOOTS- 
clear band, but is much narrower. From 329 images in 
regions covered by the SDSS, we measure a mean offset 
of 0.24 mag between the Vega-based unfiltered WOOTS 
zeropoint and the Rl zeropoint, and use it to calibrate 
the rest of the images. The scatter around this offset, 0.3 
mag, reflects both image-to-image scatter in the calibra- 
tion process (e.g., due to variations in the mix of spec- 
tral types among the lists of calibration objects), and 
systematic errors in the USNO-S photometry. The pho- 
tom etric accuracy of U SNO-.B is not well determined, 
but iMonet et a l. (2003) estimate that the scatter in the 
photometric solution over the entire catalog has an rms 
of 0.25 mag, consistent with our results. We therefore 
adopt 0.3 mag as the zeropoint error of the non-SDSS 
images. These calibration uncertainties affect the overall 
uncertainty in the SN rate, through the determination of 
our search efficiency and the visibility time (see § O ■ 

4. SN RATE CALCULATION 

Given a sample of N SNe discovered as part of a survey, 
the SN rate per unit stellar luminosity is 

N 

Ula = V T ' W 

3 

where Atj is the effective visibility time (or "control 
time"), i.e., the time during which a cluster SN Ia is 
above the detection limit of the jth image, Lbandj is 
the cluster luminosity within the search area in a chosen 
photometric band, and the summation is over all the 
survey images 1 . We explain below our measurement of 
each of these variables in the present survey. 

4.1. Visibility Time and Light Curves 

1 By "image" we will refer to the summed exposures of a field 
from a given night. 
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The effective visibility time, Atj, can be understood as 
the amount of time during which the survey is sensitive 
to SNe in a comparison between two particular images 
of a given cluster. This quantity depends on the SN 
detection efficiency in each image, the SN light curve at 
each cluster redshift, and on the frequency of searches in 
the same field. We calculate the effective visibility time 
from 

/oc 
r,*[m eff (t))dt, (2) 
-oc 

where m e ff(t) is an "effective" SN la light curve in the 
survey bandpass at the given redshift, as explained in 
§ l4.1.2l below. and ?7*[m e //(i)] is the detection probability 
as a function of SN effective magnitude. We describe 
below each step in this calculation. 

4.1.1. Detection Efficiency Estimation 

The detection efficiency of an image as a function of 
magnitude is the probability of detecting a SN-like point 
source in that image. In order to estimate the detection 
efficiency for a given field, we carry out simulations, fol- 
lowing the prescription in Gal- Yam et al. (2002). About 
200 fake SNe are added blindly to each field, with a range 
of magnitudes, and with a spatial distribution that fol- 
lows the flux of galaxies in the field. The simulated data 
undergo the same search procedure as the real data by 
means of difference image analysis (see Paper I) , and the 
number of fake SNe that are successfully detected in each 
magnitude bin is noted. The SN detection efficiency is 
usually close to 100% 1-2 mag above the detection limit 
of the image, and drops roughly linearly to zero over this 
range (Fig. [T]) . We parametrize the efficiency curve with 
the function 

!(l + e i -2 ^) , m < m . 5 
) m-m , 5 X -1 

f 1 + e s 2 J , m > mo.5, 

(3) 

where m is the Vega-based magnitude of a SN in the ef- 
fective bandpass of the survey (see below), mo.5 is the 
magnitude at which the efficiency drops to 0.5, and s 
and S2 determine the range of m over which 77 changes 
from 1 to 0.5 and from 0.5 to zero, respectively. The 
main contribution to a slow convergence to unity at mag- 
nitudes brighter than mo.5 is the difficulty of detecting 
SNe which lie close to, or are superposed on, the nucleus 
of a bright galaxy. Such cases are automatically included 
in our simulations when distributing the fake SNe in the 
images. Figure [T] shows four examples of the results of 
efficiency simulations, and the best-fit efficiency curves. 

The WOOTS database consists of 913 images, which 
were obtained under a range of observing conditions - 
atmospheric transparency, mirror reflectivity, sky back- 
ground, and telescope image quality. As a result, the de- 
tection efficiency varies significantly among the images. 
Since performing efficiency simulations in each and every 
image would be impractical, we have carried out simula- 
tions, as described above, for a subset of 30 images span- 
ning a range in observing conditions. We then searched 
for correlations between the parameters of the resulting 
efficiency curves and various parameters that character- 
ize each image. As detailed below, we found that the 



detection efficiency is primarily determined by two im- 
age parameters: the photometric zeropoint, which deter- 
mines mo.5, and the number of residuals that are left in 
each difference image, with which the upper slope param- 
eter, s, is correlated. The lower slope parameter, S2, did 
not vary significantly among the 30 simulated images, 
and was therefore set to its mean value, 0.2. 

The sensitivity of each image to SNe is determined 
by its photometric zeropoint, which combines the effects 
of detector quantum efficiency, variable mirror reflectiv- 
ity, and variable atmospheric transparency. Specifically, 
the zeropoint, ZP, is an indicator of the limiting magni- 
tude of the image, and therefore dictates the magnitude 
around which our detection efficiency drops rapidly. Fig- 
ure [2] shows the relation that we have found between ZP 
and mo.5, "*o.5 = ZP — 8.25, with a scatter of 0.42 mag. 

The number of residuals, Nn esi was defined as the 
number of objects detected by SExtractor (Bertin & 
Arnouts 1996) in each difference image, with a detec- 
tion threshold of 5er above the background. A large 
number of subtraction residuals in an image can be due 
to several reasons. Most often, poor subtraction occurs 
when the point-spread function (PSF) in the image is 
very different from the PSF in the reference image. This 
happens when the PSF is distorted, e.g., due to imper- 
fect telescope tracking, or atmospheric refraction at high 
airmass. Imperfect PSF matching and poor image sub- 
traction will mostly affect the cores of bright galaxies, 
decreasing the chances of detecting SNe that lie close to 
galactic nuclei, and causing the efficiency curve to con- 
verge more slowly to the maximal detection efficiency. 
This behavior is parametrized by s. We find that s de- 
pends on the number of residuals through the relation 
s = 0.02N Res + 0.3, with an rms of 0.3 (see Fig. [2]). 

We found no significant dependence of detection effi- 
ciency on the remaining observational parameters (sky 
background level and seeing width) due to the fact that 
most of the images were obtained during dark time, and 
the seeing spanned a limited range of 3 ± 0.4" FWHM. 
Based on the above relations, we can obtain efficiency 
curves for all images in the survey based on their mea- 
sured observational parameters, ZP and Nn es . 

4.1.2. Light Curves 

The peak magnitudes of SNe la exhibit an intrinsic rms 
scatter of 0.2-0.3 mag, and are correlated with the shape 
of the light curve through a stretch relation (Phillips 
1993; see Leibundgut 2001 for a review). Since more 
luminous SNe tend to rise and decline more slowly than 
less luminous SNe, their overall visibility time is longer 
than that of dim SNe. Sullivan et al. (2006b) study the 
dependence of SN la properties on host-galaxy type, and 
show that SNe la in elliptical galaxies tend to be dim- 
mer (with a smaller stretch factor). In their Fig. 11, 
they present the distribution of stretch factors accord- 
ing to galaxy type, at low and high redshifts. The bulk 
of the stellar light and mass in clusters is contributed 
by early-type galaxies. Indeed, among the six cluster 
SNe in our sample, the five with hosts were in early-type 
galaxies (one of the SNe was a hostless intergalactic clus- 
ter SN - see Gal- Yam et al. 2003). Assuming that all 
the SNe la in the clusters in our sample occur in ellipti- 
cal galaxies, we use the distribution of stretch factors of 
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Fig. 1. — Examples of point-source detection efficiency curves, as a function of WOOTS-clear Vega-based magnitude, for four survey 
images. Cluster names and observation dates are as marked. Circles mark the fraction of detected fake SNe, with error bars based on a 
binomial distribution. Solid curves are the best-fit efficiency curves (Eq.rjJ, and the dashed curves are the efficiency curves that correspond 
to the observational parameters of each image and the relations in Fig. 



SNe with elliptical hosts to assemble a dataset of light 
curves that will serve as templates. We base these light 
curves on a rest- frame, non-stretched , -B-band template 
light curve from iNugent et al.l (|2002f ) , and transform it 
to different stretched light curves using the stretch rela- 
tio n, M s = M s= i — a(s — l),t s = t s =\ x a as described 
by iPerlmutter et ail (|1999f h For consistency, we use the 
same non-stretched peak i?-band magnitude and a as 
Sullivan et al. (2006a), M B = -19.25 (for h=0.7) and 
a = 1.47, which are based on the results of Knop et al. 
(2003). We assume an uncertainty of 0.15 mag in Mb, 
from the dispersion in peak magnitudes of local SN light 
curves after stretch correction (Guy et al. 2005). This 
uncertainty is taken into account in our error budget (see 

§13- 

Since WOOTS observations were unfiltered, a trans- 
formation from the rest-frame .B-band light curve to an 
observed WOOTS-clear light curve at the cluster red- 
shift is necessary. For each cluster, we calculate a set 
of stretched SN light curves from multi-epoch SN la 
spectra, using synthetic photometry, as follows. We 
start with a set of m ulti-epoch spectral templates from 
INugent et al.l (|2002f) . For each combination of stretch 
and redshift, we normalize the spectra to fit the B-band 
rest-frame photometry, and redshift them according to 
the cluster redshift. We then multiply the spectra by 
the WOOTS-clear bandpass, to obtain an unfiltered light 
curve. The result is a set of representative WOOTS-clear 
light curves for every cluster, each light curve with a dif- 



ferent stretch. Figure[3]shows an example of the template 
light curves for a cluster at z = 0.131. When calculat- 
ing the visibility time for a particular image, we draw 
for it a stretch factor (with its corresponding, properly 
normalized, light curve) from the Sullivan et al. (2006) 
distribution of stretch factors. 

4.1.3. Detection Probability 

The detection efficiency function described in § 14.1.11 
is defined for the magnitude of a point source in a dif- 
ference image. This means that the interesting quantity 
in the SN light curve is not its magnitude in a single im- 
age, but the magnitude corresponding to the flux differ- 
ence A/ that results from subtracting two epochs from 
each other. This effective light curve depends on the 
time that elapsed between the two epochs that are being 
compared, 

m eff (t)=ZP- 2.5 log(A/) = (4) 
-2.51og (l0~ a4m(t) - io-°- 4m [*+(*i-*o)]^ ; 

where t is the age of the SN, and to and t\ are the times 
of the template observation and the new observation, re- 
spectively. 

The detection probability function, ri[m e ff(t)], de- 
scribes the probability of detecting a SN which occurred 
t days ago, in a comparison between two epochs. 

Since the searches in a given field were not uniformly 
distributed in time, with some observations conducted 
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Fig. 2. — Correlations between efficiency function parameters and image properties. Left: Magnitude at half the maximum efficiency 
(mo. 5) vs. image zeropoint (ZP), both in WOOTS-clear Vega magnitudes. Right: Exponential decline rate (s) vs. number of residuals in 
the difference image (Nn ea ). 

considering more than two consecutive observations, i.e., 
either 77* = 7/1 or rf = are used in equation O 

4.2. Cluster Stellar Luminosity 

4.3. Aperture Luminosity Measurement 

SN rates are often measured relative to the stellar lu- 
minosity within the search area, in a particular band. 
The cluster luminosities in our sample cannot be mea- 
sured easily from WOOTS data, since they lack color 
and spectral information from which cluster membership 
could be determined. Instead, we have used the data 
from the SDSS DR4 to measure cluster luminosities for 
72 clusters, i.e., about one-half of the sample. Unlike tra- 
ditional cluster stellar luminosity measurements, we do 
not identify the individual cluster member galaxies from 
the data. Instead, we have based our measurements on 
a method similar to aperture photometry, as follows. 
The net cluster-galaxy flux is 




-20 20 40 60 

Days since B-band maximum 

Fig. 3. — Examples of template SN la light curves with a range 
of stretch factors, calculated in the WOOTS-clear bandpass. In 
this example, the light curves are for the cluster Abell 1920, at 
z = 0.131. The representative stretch values for SNe la in low-2 
elliptical galaxies are from Sullivan et al. (2006). 

close in time to each other, the same SN can, in princi- 
ple, be discovered in more than one image. For exam- 
ple, when two closely spaced observations were compared 
to the same template, a SN that was discovered in the 
first observation (at time ti) would simply be rediscov- 
ered in the second observation (at time £2)- To avoid 
counting such objects twice in our calculation of the vis- 
ibility time, we assign the second observation a reduced 
detection efficiency, which describes the probability of 
detecting only a SN that was not detected in the previ- 
ous observation. We do this by shifting r](t) by t 2 — t\ 
days, and subtracting from it the combined detection 
probability (the probability that the SN is discovered in 
both observations): 77J = 772 — (772 x Vi)- I n principle, 
this operation can be repeated for any number of consec- 
utive observations, e.g., for three consecutive observa- 
tions, 773 will be the probability to discover a SN in the 
third observation (773), minus the probabilities of discov- 
ering it both in the first observation and in the third one, 
the second and the third, and in all three observations: 
V3 =7 l3- (V3 x 772) - (773 x 771) - (773 x 772 x 77!). In practice, 
the temporal distribution of the images does not require 



cluster 



bg 
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where k is the K-correction factor, and the summation 
is over all galaxies within a cluster radius R c that sat- 
isfy the criteria that will be described below. The total 
Galactic extinction-corrected flux 2 of galaxy i is denoted 
by ft, and f b is the average "background" -galaxy 3 flux 
per unit area. The net cluster-galaxy flux is translated to 
luminosity based on the cluster redshift and the appro- 
priate K-correction, and corrected for incompleteness, as 
described below. In what follows, we describe the deriva- 
tion of cluster stellar luminosities in the SDSS r band. 
The luminosities in three other SDSS photometric bands, 
g, i, and z, were measured in the same manner. 

To lessen the contamination of our measurements by 
foreground galaxies, we ignored, both within and out- 
side our apertures, galaxies brighter than the brightest- 
cluster galaxy (BCG) in each field. The BCG was identi- 
fied from the SDSS catalog by its magnitude, color, and 

2 Fluxes are based on the "modelmag" SDSS magnitudes. 

3 Non-cluster galaxies are, of course, both foreground and back- 
ground, but we retain the term in analogy to aperture photometry. 
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where available, by redshift. We also ignored galaxies 
fainter than the SDSS 95% completeness limit (see Ta- 
ble [2J. 



TABLE 2 
Adopted Photometric Parameters 



from the faint end of a Schechter (1976) luminosity func 
tion: 

J™HL)dL 



C = 





(1) 


(2) 


(3) 


(4) 


band 


m lim 


A^band 


M* 


a 




[mag] 


[mag] 


[mag] 




9 


22.2 


5.12 


-22.01 ± 0.26 


-1.00 ±0.06 


r 


22.2 


4.64 


-22.21 ± 0.05 


-0.85 ± 0.03 


i 


21.3 


4.53 


-22.31 ± 0.08 


-0.70 ±0.05 


z 


20.5 


4.51 


-22.36 ± 0.06 


-0.58 ±0.04 



Note. — (1) SDSS 95% completeness magnitude; Adelman- 
McCarthy et al. (2006). 

(2) Solar absolute magnitudes; Blanton et al. (2006). 

(3) , (4) Schechter luminosity function parameters for clusters; Goto 
et al. (2002). 



Leir & van den Bergh (1977) estimated cluster radii, 
RlvABi f° r 1889 Abell clusters, using the red plates of 
the Palomar Sky Survey, by superposing a grid of circles 
of various radii on each cluster, and selecting by eye the 
radius that encompasses most uniformly all visible clus- 
ter members. We have examined their estimate for 12 of 
our clusters, by measuring the net cluster-galaxy flux as 
described above, as a function of R c . We find that the 
enclosed flux profile generally tends to a constant value 
(to within errors) at radii R c « RlvclBi an d we there- 
fore adopt RLvdB as the cluster radius of each cluster 
in our sample. In §3 we investigate the effect of vary- 
ing the choice of R c between 0.5Rl V cIB and l.5RLvdB, 
according to a normal distribution centered on RlvcIB 
with a — 0.2RL V dB- We find that this uncertainty in 
the cluster radius results in a ~ 1% error in the stellar 
luminosity. 

The average background-galaxy flux per unit area, f bg , 
was estimated from the flux per unit area of galaxies that 
satisfy the above criteria, in an annulus of area of ~ 0.5 
deg 2 centered on the cluster, with inner radius of 0.4°. 
The inner radius is chosen to be large enough to exclude 
cluster member galaxies, but small enough to be repre- 
sentative of the large-scale structure in the vicinity of the 
cluster. In order to check whether the measured cluster 
luminosity is sensitive to the choice of background an- 
nulus inner radius, we calculated the luminosity of each 
cluster using a range of inner radii, between 0.35° and 
0.45°. We experimented with both using a fixed angular 
radius for all of the clusters, and using a buffer that is 
proportional to the cluster radius. The luminosities that 
result from the various methods of background estima- 
tion have an rms of 8%. We take this uncertainty into 
account in our error budget (see fj5]). 

To convert the observed magnitudes to rest-frame mag- 
nitudes, a K-correction factor was calculated for each 
cluster assuming that all the cluster light is emitted by 
elliptical galaxies at the cluster redshift. The net cluster 
flux was then translated to luminosity, according to the 
cluster redshift and the adopted cosmology. The faint 
cut on the galaxy magnitudes means that the summed 
luminosity includes only galaxies brighter than some lim- 
iting luminosity, which depends on the cluster redshift. 
We correct for this incompleteness by multiplying the 
summed luminosity by the fraction of light that comes 



roc 



(mn m ,z) 



§{L)dL 



(6) 



where $(L)dL = (L / L*) a ex.p(—L / L*)d(L / L*) . We 
adopt a = -0.85 ± 0.03 and M* = -22.21 ± 0.05 as the 
mean values for the r-band luminosity function parame- 
ters in clusters (Goto et al. 2002; see Table[5]). The inte- 
grated r-band luminosity correction factors in our sample 
are in the range C-l = 0.4±° °? x 10~ 3 to 8.3tj^ x 10~ 3 , 
for z = 0.06 and z = 0.2, respectively, i.e., always quite 
small (see Table [3]). We note that Goto et al. (2002) 
excluded the BCGs from the analysis when fitting the 
data to the Schechter function. We therefore applied the 
correction above after subtracting the luminosity of the 
BCG from the summed luminosity, and added the BCG 
luminosity to the corrected cluster luminosity. 

Several recent studies have argued that BCG galaxy 
luminosities are underestimated in SDSS, because the 
default sky subtraction algorithm removes the outer, low 
surface-brightness, flux from these galaxies (Graham et 
al. 2005; Bernardi et al. 2005; Lauer et al. 2006; 
Desroches et al. 2006; von der Linden et al. 2006). L.- 
B. Desroches and C.-P. Ma (private communication) find 
that, for a subsample of the Miller et al. (2005) C4 sam- 
ple, improved sky subtraction increases the luminosity 
by a mean factor 1.3 over the total deVaucoulers SDSS 
luminosity. For five of our BGCs which are in common 
with their subsample they find similar correction factors. 
It is clear that BCG luminosity measurements from the 
SDSS are indeed sensitive to the sky-level subtraction al- 
gorithm. On the other hand, it is not certain that the 
new corrected measurements are in fact superior, e.g., 
they might include in the BCG measurement some light 
from neighboring galaxies. We therefore choose to in- 
crease the BCG luminosities by 15%, with an additional 
systematic uncertainty of ±15%. Since the BCG lumi- 
nosity typically constitutes about 6% of the cluster lumi- 
nosity, this correction lowers our derived SN rate by 1%, 
and increases the systematic error in the rate by ±1%. 

4.3.1. Comparison to Other Luminosity Measurements 

The reliability of our luminosity measurement method 
can be tested by comparing our results to independent 
measurements of the same clusters. The comparison can 
be particularly straightforward, if the compared lumi- 
nositie s were deriv e d from the same data, namely, the 
SDSS. iMiller et all (|2005l ) describe their compilation of 
a cluster catalog, containing 748 clusters selected spec- 
troscopically from the SDSS database, 12 of which are 
in our samp l e and are fully covered by the SDSS DR4. 
IMiller et alj (|2005l ) calculate the r-band luminosities of 
the clusters as the summed luminosities of all galaxies 
within 2.1 Mpc (for our choice of Hq) from the center 
of the cluster, within 4er in redshift space, and more lu- 
minous than 1.4 x 10 10 L Q (corresponding to a 17.8 mag 
galaxy at z = 0.11). In order to compare our results to 
theirs, we have calculated the aperture-based luminosi- 
ties for the 12 clusters that appear in both samples, with 
the same radius and limiting luminosity criteria. Exclud- 
ing Abell 1539, for which the luminosity measurement of 
Miller et al. was corrupted by a satellite trail in the SDSS 
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TABLE 3 

Measured Cluster Luminosities 



Abcll 
Number 


z 




9 


Lu 




osity 


[10 ij! 

i 


1-0] 
z 


L 

C 


uminosity 
orrcction (r) 


Mass 

[10 12 A* Q ] 


125 





19 


2 


51 


2 


51 


3.10 


4.50 
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1 
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1 
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1 


55 


2 


06 


2.32 
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1 


0061 
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94 
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1 


0022 
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2 


81 


3 


30 


3.77 


4.77 


1 
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4.8 ±1.1 
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2 


27 


3 


18 


3.46 
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1 


0054 
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1227 





11 


2 


45 


3 


19 


4.50 
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17 
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11 
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32 
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17 
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46 


5 


39 
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16.1 ±3.8 
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4.83 
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1 
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r-band image, we find that the mean difference between 
the two measurements of every cluster is ~ 20%. 

4.3.2. Measured Luminosities 

Table [3] lists our luminosity measurements for the 72 
clusters in our sample that are covered by the SDSS DR4 
(the "SDSS clusters"). Due to the WOOTS selection 
criteria (Paper I), the luminosities span a small range, 
with 68% of the clusters having L r — 2.3+J;g x 10 12 £ riQ 
(Fig. [4|. Within this range, there is no clear trend with 
other cluster parameters. We therefore assign to the clus- 
ters for which we do not have SDSS data (the "non-SDSS 
clusters" ) random luminosities drawn from the luminos- 
ity distribution of the SDSS clusters. We estimate the 
uncertainties in the SN rate due to this random luminos- 
ity assignment in § [SJ 

The actual luminosity that we require in Eq. [4] is not 
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Fig. 4. — Distribution of stellar luminosities of the WOOTS clus- 
ters that are in the SDSS DR4 (the "SDSS clusters"), as measured 
with our photometric aperture technique. 

the total cluster luminosity, but the cluster luminosity 
included within the search area. For every WOOTS im- 
age we define the effective search area as the overlap area 
between the image and its reference image. The reference 
image was assigned a zero search area, since the survey 
was sensitive only to SNe that are brighter in the later 
epoch (see Paper I). For the SDSS clusters, we calculate 
the luminosity, as described above, contributed by galax- 
ies that are inside the effective area of an image pair. If 
the full cluster area within R c is covered by the image 
(usually clusters with R c < 6'), we adopt the calculated 
cluster luminosity instead. 

To estimate the luminosity inside the search area of the 
non-SDSS clusters, we need to scale down the cluster lu- 
minosity according to the luminosity profile and the effec- 
tive search area. Figure [5] shows, for the SDSS clusters, 
the mean fraction of cluster luminosity, L(< r)/L(< R c ), 
measured within an aperture of radius r, vs. the radius 
normalized by the cluster radius, r/R c . The solid line 
is a second-order polynomial fit to the data, L r /LiR c — 
-0.07+ lMr/R e - 0.56(r/i?c) 2 . We use this relation for 
all the non-SDSS clusters with R c > 6', adopting for r 
the radius of a circle with an area equal to the effective 
search area of the image. Figure [6] examines the relia- 
bility of this approximation. For each SDSS cluster, we 
estimate the luminosity within a circular aperture with 
an area of 12' x 12' (the dimensions of the CCD used in 
WOOTS), using the above relation, and plot it against 
the luminosity that is measured directly in the effective 
search area. The difference between the estimated lumi- 
nosities and the measured ones has an rms deviation of 
16%, which we adopt as the uncertainty introduced by 
this approximation. 

5. RESULTS AND ERROR ESTIMATION 

Combining all the elements of the calculation, as de- 
scribed above, we can derive the cluster SN rate and 
its uncertainty. The dominant source of error in the 
rate is the Poisson statistics of the (small) number of 
SNe, N = 6±|:H (68% confidence limit). To determine 
the propagation of the uncertainties in all the param- 
eters and measurements that enter the rate derivation 
to the overall systematic uncertainty, we have conducted 
a Monte Carlo experiment, in which we calculate the 
SN rate many times, with the parameters drawn at ran- 
dom each time from their respective distributions. The 
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Fig. 5. — Composite enclosed cluster-galaxy light profile, based 
on all SDSS clusters in the sample. Error bars represent the mea- 
sured rms over the sample. The solid line is the indicated second- 
order polynomial fit to the data. 
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Measured luminosity [10 L^] 

Fig. 6. — Comparison between the luminosity within the 
WOOTS search area, estimated using the relation in Fig. \E\ and 
the actual luminosity measured in the image area, for the SDSS 
clusters. The solid line represents perfect agreement. 

distributions are either Gaussian with an rms deviation 
set to equal the parameter error, or a specific measured 
distribution (e.g., the distribution of light-curve stretch 
factors). The distributions of SN rates from the Monte 
Carlo simulation permit both a reliable error estimation 
and a determination of the most probable rate. Be- 
cause of the nonlinear dependence of the rate on some 
of the parameters entering the calculation (e.g., stretch 
factor and luminosity), the most probable SN rate will 
not necessarily be the SN rate obtained from the com- 
bined most-probable values of all the input parameters. 
We will take as the most probable SN rate the peak of 
the Monte Carlo distribution, with a systematic error 
derived directly from this distribution, and a statistical 
error that originates from the Poisson confidence limits 
on the number of events. 

We have also examined the sensitivity of the SN rate to 
the uncertainties in the individual parameters, by turn- 
ing on the Monte Carlo simulation for each one sepa- 
rately. Table [4] summarizes the mean values and the un- 
certainty in the SN rate, due to each of the parameters. 
Luminosity error: The overall luminosity error is 
about 5%, and is dominated by the process of randomly 
drawing the non-SDSS cluster luminosities from the lu- 
minosity distribution of the SDSS clusters. Varying the 
cluster radius has a much smaller effect, of ~ 1%. Sim- 



ilarly, the uncertainty due to the luminosity gradient of 
the clusters introduces an SN rate error smaller than 
0.2%. Varying the inner radius of the background an- 
nulus results in an error of less than 0.5%. 
Visibility time error: The uncertainties in the visi- 
bility time measurement have a less intuitive effect on 
the final rate distribution. The visibility time is most 
sensitive to the shape and parameters of the efficiency 
function, mainly mo. 5, which is the magnitude at which 
the efficiency drops to half its maximum value. We find 
that drawing the efficiency parameters from normal dis- 
tributions centered on their best estimated or mean val- 
ues results not only in a dispersion in SN rates, but also 
in an overall decrease in the rate. The reason for this 
behavior is that the dependence of the rate on the effi- 
ciency is nonlinear. A slightly higher efficiency (caused 
by higher mo. 5 or lower s) will result in a decrease in the 
SN rate, and vice versa, but for a similar absolute change 
in the efficiency parameters, the decrease in the rate is 
larger than the increase in the rate. The same applies to 
the distribution of SN light-curve stretch factors. Draw- 
ing light curves from a stretch-factor distribution, even if 
it is symmetric about the mean value, causes an overall 
decrease in the mean SN rate. In our case, all of these 
effects are, of course, negligible compared to the statisti- 
cal error. The total systematic error, which in our case 
is ~ 6%, will become comparable to the statistical error 
in future surveys with the detection of several hundred 
SNe. 

We note that some parameters were not varied in our 
simulation, since their errors arc significantly smaller. (1) 
The uncertainties in the Schechter luminosity function 
parameters cause a symmetrical change of 0.16% in the 
luminosity, in the most extreme case, and were therefore 
ignored. (2) In the calculation of cluster luminosities, 
we assumed that all of the cluster light is emitted by el- 
liptical galaxies. However, a small fraction of the light 
does come from spiral galaxies, which have a different 
K-correction - it is 0.16 mag smaller at z = 0.06, and 
0.26 mag larger at z = 0.2. To check the effect of this 
assumption, we calculated the SN rate for the unrealistic 
case, in which half of the cluster galaxies are Sc galaxies. 
In the 0.06 < z < 0.2 range, individual cluster luminosi- 
ties increase (at z < 0.115) or decrease (at z > 0.115) by 
no more than 10%, relative to the all-elliptical case. The 
actual distribution of cluster redshifts leads to a ~ 4% 
increase in the SN rate. Thus, even in this extreme sce- 
nario, the influence on the SN rate is smaller than that 
of the other sources of systematic error. 

The resulting SN rates, as determined from the po- 
sitions of the peaks of the Monte Carlo distributions, 
and normalized by stellar luminosity in the g,r,i, and z 
bands are given in Table [5l The rates are given in units 
of SNuband, defined as SNe (century 10 10 Lband,o) -1 - We 
also express our result in the Johnson B band used tra- 
ditionally for SN rates, in units of SNus, by assum- 
ing B — r = 1.078 mag for elliptical galaxies (as cal- 
culated from the template elliptical spectrum of Kinney 
et al. 1996), M rQ = 4.64 (Blanton et al., 2006), and 
M fl)0 = 5.47 (Allen 1976). 

6. COMPARISON TO OTHER SN RATE MEASUREMENTS 

Our derived cluster SN rate at 0.06 < z < 0.19 can 
be compared to recent measurements of SN rates both 
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TABLE 4 

Sensitivity of the SN Rate to the Uncertainty in Individual Parameters and Measurements. 





Parameter distribution 
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S 1 IKS Hust or 


0.325 


0.016 




luminosity distribution 






Cluster radius (i? c ) 


Normal with a = 20% 


0.327 


0.004 


L(in image) oc L(< r) 


Normal with a = 0.16 


0.324 


0.002 


Background annulus 


Measured 1 ; a = 0.07 


0.3243 


0.0007 


Overall uncertainty in L 




0.327 


0.017 


ZP (from USNO-B / SDSS) 


Normal with a = 0.3/0.04 


0.316 


0.004 


m0.5Ds.ZP 


Normal with a = 0.4 


0.301 


0.007 


svs.N Res 


Normal with a = 0.3 


0.321 


0.002 


LC stretch factor 


From Sullivan et al. (2006b) 


0.318 


0.004 


Overall uncertainty in visibility time 




0.287 


0.009 


Overall uncertainty from all parameters 




0.290 


0.017 



By calculating the luminosity for different buffer selections; see § 14.21 



TABLE 5 
SN Ia Rates 



Value 1 



Units^ Reference 



Environment Redshift 



cluster 



0.06 < z < 0.2 0.351 



±0.020 
±0.018 
±0.014 
±0.010 



0.139 

0.288± '^2 

0- 2 29to.o1I 
186+° 111 

U.iOD_ f) Q74 

0.36to^4 ±0.02 

n+0.59 
-0.25 

O.BOtfcg 
0.16 ±0.05 



SNu 9 
SNn r 
SNui 
SNu z 



This work 



cluster 
cluster 
cluster 
E/SO 



0.06 < z < 0.2 
0.18 < z < 0.37 
0.83 < z < 1.27 
local 



0.39" 1 



0.098loo39 ± °- 009 
0.038ir « 
0.053 ± 0.011 



SNu B This work 

SNu s Gal- Yam et al. (2002) 

SNu s Gal- Yam et al. (2002) 

SNu B Capellaro et al (1999) 



cluster 

E/SO 

E/SO 



0.06 < z < 0.2 
local 

0.2 < z < 0.75 



SNuM This work 

SNuM Mannucci et al. (2005) 

SNuM Sullivan et al. (2006) 



in clusters and in the field. In the field, we compare 
the rate derived in this work to recent measurements 
of the SN rate in elliptical galaxies per unit stellar lu- 
minosity and per unit mass. Cappellaro et al. (1999) 
measured a local E/SO SN Ia rate of 0.16 ± 0.05 SNu B , 
which is consistent with the rate we obtain in this lumi- 
nosity band, 0.36^^4 ±0.02 SNub. To convert our mea- 
surement to SNuM [SNuM = SN (century 1O 1O M )- 1 ], 
we follow Mannucci et al. (2005), who used the mass- 
to-light ratio derived by Bell et al. (2001) to con- 
vert the ii-band luminosities and B — K colors of the 
galaxies in their sample to stellar mass. We derive 
the r — i color of each SDSS cluster from the ratio 
between its measured r-band and z-band luminosities, 
r - i = -2.51og 10 (L r /i j; ) - M i:Q + M r>Q . The stellar 
mass of each cluster is then estimated from the color- 
dependent stellar mass-to-light ratio derived by Bell et 
al. (2003), log lo (M /L z , Q ) = -0.052 + 0.923(r - i) (see 
Manucci et al. 2005 for a discussion of the validity of 
this ratio for our purpose). Finally, we define the stel- 
lar mass-to-light ratio of our cluster sample as the total 
mass divided by the total z-band luminosity of the entire 
sample, M/L = 1.89 ± 0.44M Q /L 2iQ , and use it to con- 
vert the SN rate from SNu z to SNuM. The resulting SN 
rate per unit stellar mass is 0.098^^^ ± 0.009 SNuM. 
Mannucci et al. (2005) used the SN sample of Cappel- 
laro et al. (1999) to measure the SN rate per unit mass 
as a function of host-galaxy morphological type. They 
found a rate of 0.0383q'q^2 (converted from h = 0.75 to 
h = 0.70) SNuM for low-redshift (2 < 0.02) E/SO galax- 



ies, lower than, but consistent with our result for cluster 
galaxies at a somewhat higher redshift. We note that 
some of the E/SO galaxies monitored by Cappellaro et al. 
are members of nearby galaxy clusters, and therefore the 
Mannucci et al. (2005) local rate includes both cluster 
and field early- type galaxies. A discussion of local rates 
separated by environment will be presented elsewhere. 

Sullivan et al. (2006) measured the SN Ia rate at 
0.2 < z < 0.75, from 124 SNe discovered by the Su- 
pernova Legacy Survey (SNLS), and studied the SN 
host properties. Following Mannucci et al. (2005) 
and Scannapieco et al. (2006), they separated the rate 
into two components, one proportional to stellar mass, 
and the other proportional to the star formation rate - 
SNR Ia (t) = AM(t) + BM(t). Their best-fit values are 
A = 0.053 ± 0.011 SNuM, B = (3.9 ± 0.7) x 10~ 4 SNe 
yr _1 (M Q yr -1 ) -1 . Since the stellar mass in their sam- 
ple is dominated by the early-type galaxies, and these 
galaxies have virtually no star formation, the parameter 
A is essentially the SN rate in E/SO galaxies. Compar- 
ison of the E/SO SN rates measured by Mannucci et al. 
(2005) and Sullivan et al. (2006) indicates that the rate 
may be constant with redshift out to z « 0.5, although 
an increase by a factor ~ 2 is also consistent within the 
errors. Our measurement in clusters at an intermediate 
redshift range is consistent with these two field elliptical 
rate measurements. 

In galaxy clusters, Gal- Yam et al. (2002) measured 
cluster SN Ia rates per unit stellar B-band luminosity of 
°- 39 to.25 SNu s at < z >= °- 25 and O-SOtoii SNu s at < 



SN la Rate in Low-Redshift Galaxy Clusters 



11 



z >= 0.9 (converting to our adopted cosmology). These 
estimates were based on the detection of one cluster SN in 
the low-redshift bin, and one or two in the high-redshift 
bin, and consequently, the errors are large. Comparing to 
our measurement of 0.36^° '^±0.02 SNub, there is a hint 
for a slow rise in SN rate with redshift. However, the data 
are equally consistent, given the large error bars, with a 
constant rate. More accurate SN rate measurements at 
high z, which are in progress, will elucidate this question. 

7. SUMMARY 

We have analyzed the data from the WOOTS cluster 
SN survey (Paper I) to derive the SN la rate in 0.06 < 
z < 0.19 clusters, normalized by stellar luminosity and 
by stellar mass (Table [5]). In addition to the sample of six 
cluster SNe la discovered by the survey, our measurement 
required the determination of two variables - the survey 
visibility time, and the cluster stellar luminosities, which 
wc determined using the survey data, assisted by SDSS 
data. We have conducted Monte Carlo simulations that 
mimic the inhomogeneity in the parameters that enter 
the rate calculation, and quantify the dependence of the 
systematic errors on the uncertainty in these parameters. 
The resulting systematic errors are about an order of 
magnitude smaller than the Poisson errors due to the 
small number of SNe that were detected. 

We find that the SN rate is similar to the rate found 



in field elliptical galaxies, both locally and out to z = 
0.75. A comparison to cluster SN rates at higher redshifts 
is similarly consistent with an unchanging rate, but the 
large uncertainties in current high-z cluster SN rates can 
also accommodate a decline in the rate with time. We are 
in the process of obtaining cluster SN rate measurements 
at high redshift. The emerging dependence of SN rate on 
cosmic time and on galaxy environment should provide 
valuable new insights for astrophysics and cosmology. 
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